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Analytic versions of QCD are those whose coupling aa{Q^) does not have the unphysical Landau 
singularities on the space-like axis (— g^ — Q^ > 0). The coupling is analytic in the entire complex 
plane except the time-like axis {Q^ < 0). Such couplings are thus suitable for application of per- 
turbative methods down to energies of order GeV. We present a short review of the activity in the 
area which started with a seminal paper of Shirkov and Solovtsov ten years ago. Several models 
for analytic QCD coupling are presented. Strengths and weaknesses of some of these models are 
pointed out. Further, for such analytic couplings, constructions of the corresponding higher order 
analytic couplings (the analogs of the higher powers of the perturbative coupling) are outlined, and 
an approach based on the renormalization group considerations is singled out. Methods of evalu- 
ation of the leading-twist part of space-like observables in such analytic frameworks are described. 
Such methods are applicable also to the inclusive time-like observables. Two analytic models are 
outlined which respect the ITEP Operator Product Expansion philosophy, and thus allow for an 
evaluation of higher-twist contributions to observables. 

PACS numbers: 12.38.Cy, 12.38.Aw,12.40.Vv 
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I. INTRODUCTION 

Perturbative QCD calculations involve coupling 
a{Q^) = as{Q'^)/TT which has Landau singularities (poles, 
cuts) on the space-like semiaxis < Q^ < A^ {q^ = 
— Q^). These lead to Landau singularities for the evalu- 
ated space-like observables ViCp) at low Q^ ~ A^. The 
existence of such singularities is in contradiction with the 
general principles of the local quantum field theories \M, ■ 
Further, lattice simulations [2] confirm that such singu- 
larities are not present in a{Q^). 

An analytized coupling Ai{Q'^), which agrees with the 
perturbative a{Q^) at Q^ ^ oo and is analytic in the Eu- 
clidean part of the Q^-plane {Q^ ^ C, Q"^ ^ 0), addresses 
this problem, and has been constructed by Shirkov and 
Solovtsov about ten years ago 0]. 

Several other analytic QCD (anQCD) models for 
Ai{Q'^) can be constructed, possibly satisfying certain 
additional constraints at low and/or at high Q^. 

Another problem is the analytization of higher power 
terms a" i— >• An in the truncated perturbation series 
(TPS) for 'D{Q^). Also here, several possibilities appear. 

Application of the Operator Product Expansion (OPE) 
approach, in the ITEP sense, to inclusive space-like ob- 
servables appears to make sense only in a restricted class 
of such anQCD models. 

This is a short and incomplete review of the activity 
in the area; relatively large space is given to the work of 
the review's authors. For an earlier and more extensive 
review, see e. g. Ref. [J|. 



Section |TT] contains general aspects of analytization of 
the Euclidean coupling a(Q^) i-^ Ai{Q^), and the def- 
inition of the time- like (Minkowskian) coupling ^i(s). 
Further, in Sec.|TT]we review the minimal analytization 
(MA) procedure developed by Shirkov and Solovtsov [31] , 
and a variant thereof developed by Nesterenko @- In 
Sec. mil we present various approaches of going beyond 
the MA procedure, i.e., various models for 2li(s), and 
thus for AiiQ^) [i,0,S,i9, l_0,li]- In Sec.lvl analytiza- 
tion procedures for the higher powers a"((5^) >-^ AniQ"^) 
in MA model are presented [l2, llj, UM i ^^'^ ^^ alterna- 
tive approach which is applicable to any model of ana- 
lytic Ai{Q'^) [13, Ell is presented. In Sec. |Vl an analy- 
tization of noninteger powers a'^{Q^) is outlined [ISJ. In 
Sec. Ell methods of evaluations of space- like and of inclu- 
sive time-like observables in models with analytic Ai (Q^) 
are described, and some numerical results are presented 
for semihadronic r decay rate ratio r^, Adler function 
dv(Q^) and Bjorken polarized sum rule (BjPSR) db[Q^) 
m> El, US, [llHJIl- In Sec. Eni two sets of models 
are presented [13, 13 whose analytic couphngs Ai{Q'^) 
preserve the OPE-ITEP philosophy, i.e., at high Q^ they 
fulfi ll: \A i[Q^) - a(Q2)| < [A?/Q-^f for any k e N. Sec- 
tion IVIIII contains a summary of the presented themes. 



II. ANALYTIZATION a(g2) H^yii(g2) 

In perturbative QCD (pQCD), the beta function is 
written as a truncated perturbation series (TPS) of cou- 
pling a. Therefore, the renormalization group equation 
(RGE) for aiQ"^) has the form 
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da{\nQ^-P2,--.) 
91nQ2 



^/3j_2a^(lnQ2;/32, 

J=2 



(1) 



The first two coefficients [/3o = (1/4)(11 - 2n//3), Pi = 
(1/16)(102 — 38n^/3)] are scheme-independent in mass- 
independent schemes. The other coefficients (/32,/33, . . .) 
characterize the renormahzation scheme (RSch) . The so- 
lution of perturbative RGE ([T]) can be written in the form 
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(2) 



where L = ln(Q^/A^) and Kf^^ are constants depending 
on /3/s. In MS: A = A - lO"! GeV. 

The pQCD couphng a{Q^) is nonanalytic on — oo < 

Q^ < A . Apphcation of the Cauchy theorem gives the 
dispersion relation 
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(ry -> 0), (3) 



where pf {a) is the (pQCD) discontinuity function of a 
along the cut axis in the Q^-plane: pf {a) — Ima(— cr — 
ie). The MA procedure of Shirkov and Solovtsov [3| 
removes the pQCD contribution of the unphysical cut 
< — cr < A^, keeping the discontinuity elsewhere un- 
changed ( "minimal analytization" of a) 
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In general: 
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(5) 



where pi{cr) = Im^i(— cr — ie) . Relation (O defines 
an analytic coupling in the entire Euclidean complex Q^- 
plane, i.e., excluding the time- like semiaxis —s = Q^ < 0. 
On this semi-axis, it is convenient to define the time-like 
(Minkowskian) coupling 2li(s) (Hill, [11 
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The following relations hold between ^i, 2li and pi: 



2ti s = - / —pi{a) , 
ds%iis) 
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(6) 

(7) 
(8) 
(9) 



The MA is equivalent to the minimal analytization of the 
TPS form of the /?(a) = da{Q^)/d\nQ^ function [l^ 



dAi^''^HlnQ^;P2,...) 
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where pg (cr) = Im/3(a)(— cr — ie), and 



Pia) 



-^/3,_2a^(lnQ2;^2,...). (n) 

The MA couplings AiiQ"^) and 2ti(s) are finite in the 
IR (with the value l//3o at Q^ = 0, or s = 0) and 
show strong stability under the increase of the loop- 
level ri„i — jniax — 1 (see Figs. [U [2]), and under the 
change of the renormahzation scale (RScl) and scheme 
(RSch). Another similar pQCD-approach is to analy- 
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FIG. 1: Left: one-loop MA oeCQ) = i"^i(Q^) and its one-loop 
perturbative counterpart as{Q'^) in MS, for Uf = 3 and A = A = 
0.2 and 0.4 GeV. Right; stability of the MA «e(Q) = nAi(Q^) 
under the loop-level increase. Both figures from: Shirkov and 
Solovtsov, 1997 H. 
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FIG. 2: The MA time-like and space-like couplings 2li(s^") and 
^i(Q) at 1-loop, 2-loop (3-loop) level; in MS for nf = 3 and 
A = 0.35 GeV [2li and ^i in figure are ttSIi and ttA^i in our nor- 
malization convention]. Figure from: Shirkov and Solovtsov. 2006 

tize minimally Pia)/a = d \n a{Q^) / d \n Q"^ [J, B IIH. 
This leads to an IR-divergent analytic (MA) coupling, 
A(Q') -^ (AVQ2)Qn(AVQ'))-' when Q^ -^ 0. At 
one-loop: 
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/3o(OVA2)ln(Q2/A2) 



(12) 



Also this coupling has improved stability under the loop- 
level change, and under the RScl and RSch changes (see 



Figs. [3l HI). Numerical predictions of this model, at the 
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FIG. 3: Left: one-loop MA aan(Q) = f3oAi{Q^) and its one- 
loop perturbative counterpart, as a function of Z = Q^/A^ (Fig- 
ure from: Nesterenko, 2000 [3|). Right: stability of the MA 
5an(Q) = /3o^l(Q^) under the loop- level increase, as a function 
of Z = Q2/yv2 (Figure from: Nesterenko, 2001 ^). 
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FIG. 4: One-loop time-like and space-like MA couplings &a,u(s) = 
7r2li(s) and Oan(Q^) = ''r.4i(Q^) as a function of Z = —s/h? or 
Z = Q^/A^, respectively. Figure from: Nesterenko, 2003 |21| . 

one-loop level, for various observables, were performed 
in Ref . [2l| , and they agree with the experimental results 
within the experimental uncertainties and the theoretical 
uncertainties of the one-loop approximation. 



III. BEYOND THE MA 

The idea to make the QCD coupling IR finite phe- 
nomcnologically is an old one, by the substitution 
ln((5^/A^) K^ ln[((5^-l-4mo^/A^] where nig is an effective 
gluon mass, cf. Refs. [13,0, [11]. 

On the other hand, the analytic MA, or MA, couplings 
can be modified at low energies, bringing in additional 
parameter(s) such that there is a possibility to reproduce 
better a wide set of low energy QCD experimental data. 

Among the recent proposed analytic couplings are: 

1. Synthetic coupling proposed by Alekseev [6|: 



a,y„(Q2) = a(MA)(Q2)^ 



/3o 



cA2 



rfA2 



Q2 



, (13) 



where the three new parameters c, d and gluon mass 
nig were determined by requiring a^yniQ'^) — otptiQ^) ^ 
{K^ IQ'^Y (foi' tliG convergence of the gluon condensate) 
and by the string condition V{r) '^ ar {r ^ oo) with 
a ~ 0.42^GeV^. This couphng is IR-divergent. 
2. The coupling by Sriwastawa et al. [7|: 
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(14) 



{a + z ~ ie){a + l){l + zP) 



where z — Q^ / h? and < p < 1. This formula coincides 
with Ncstcrcnko's (one-loop) MA coupling when p — 1- 
3. An IR- finite coupling proposed by Webber [§]: 
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(15) 



where z = Q^/A^ and specific values are chosen for pa- 
rameters b = 1/4, c = 4, and p = 4; ajy (0) ~ 7r/(2/3o). 

4. "Massive" MA or MA couplings AiiQ^) and 2li(s) 
proposed by Nesterenko and Papavassiliou [^]: 
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Q2 _|_ 4„j2 
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a — Am da 



(16) 

where m ^ A; and pi{a) = pi^^^\a) in the MA case. In 
this case: ^i("''(0) = Sti^'^^O) = 0. The mass to is some 
kind of threshold, and can be expected to be ^ m^^. 

5. Two specific models of IR- finite analytic coupling 
[lOl . [ll|: on the time-like axis s = —Q^ > 0, the par- 
turbative discontinuity function pi(s), or equivalently 

Sli*- '(s), was modified in the in the IR regime (s ~ A ). 
A first possibility (model 'Ml'): 

2l(^i)(s) = CfMls{s-Ml) 

+fcoe(Mo - s) + e(s - Mo)2tl"^^ (s) , 

2 2 2 2 

where c/, fcg, Cr = Af^/A , cq = Mq/A are four dimen- 
sionless parameters of the model, all ~ 1. One of them 
(/co) can be eliminated by requiring the (approximate) 
merging of Ml with MA at large Q^: 

\A,^^'\Q')-A^^''^\Q')\^{A'/Q'r. 



The Euclidean Af^^\Q^) is 



l(Ml) 
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(17) 



where the constant df is 

Another, simpler, possibihty is (model 'M2'): 

str'w = 2ir^'w+c.e(Mj-.), (18) 



M„ 



Ar\Q') = Ar'^'m+c. '_, , (w) 



{Q'+Ml) 



72 ,-r2 



where c« and Cp = Mp/A are the model parameters. 

6. Those anQCD models which respect the OPE-ITEP 
condition are presented in Sec. I VIII 



analogously for dAi/dPs, etc. In our approach, the basic 
space- like quantities are Ai{fi^) of a given anQCD model 
(e.g., MA, Ml, M2) and its logarithmic derivatives 



AnifJ-'^) 



2^ _ ^-ir-^9-Ui(M!) ,^,„ ^ 

/3J'-i(n-l)!9(ln^2)„-i ' ^^ ^,^,---), 



whose pQCD analogs are 



(23) 



°"^^ ^^/3o"-i(n-l)!a(lnM2)"-i ' ^'^ " 1' 2' • • •)• 

(24) 
At loop- level three {n^ = 3), where we include in RGE 
dll) term with jmax = 4 (thus /32), relations (PT|) are 



IV. ANALYTIZATION OF HIGHER POWERS 

a' 1-^ Ak 

In MA model, the construction is [1 [H, [H, [11 

(MSSSh: Milton, Solovtsov, Solovtsova, Shirkov): 



a'(Q') 
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(MA), ^2 
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da 
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pri'^) : (20) 



where /c = 1,2, . . .; p^'' ^(ct) = Im[a'^(— cr — ie)]; and a is 
given, e.g., by Eq. Q. In other words, "minimal analy- 
tization" (MA) is applied to each power a''. 
As a consequence, in MA we have [19] 

a <(MA)/ 2\ 

o m /i^ 

02 /I (MA) / 2\ 

^ -^1 y^^ ) _ 9fl2/l(MA) .^ ^ ^(MA) 

gn^ 2)2 - ^Po-^a +5/:;oPiA +•••. 

etc. This is so because a'^, and consequently pj^ (c), 
fulfill analogous RGE's. 

The approach (PD)) of constructing ^^'s (fc > 2) can 
be applied to a specific model only (MA). In other an- 
QCD models (i.e., for other AiiQ"^)), the discontinuity 
functioiis Pk (k > 2) are not known. We present an ap- 
proach 10, 11] that is applicable to any anQCD model, 
and reduces to the above approach in the MA model. We 
proposed to maintain the scale (RScl) evolution of these 
(truncated) relations for any version of anQCD 



dAi{p^;P2,...) 

9Ml(//2;/32,...) 



a(ln^2)2 



— —fioAl — ■ • • — f3n^-2An^ , 

= 2/?2^3 + 5/3o/3i^4 + --- + 4!l^„„, 

(21) 

etc. Eqs. dH]) define the couphngs A(Q^) (fc > 2). Fur- 
ther, the evolution under the scheme (RSch) changes will 
also be maintained as in the MA case (and in pQCD): 



^Al{^i'■,f32,■■■) 

d(32 



(22) 



^2(m') = A(/i') + ^^3(M'), Asi^i') = Asi^i'), (25) 



/3o 



implying 



^2(m') = A2{^l')-^A^i^i^), Asifi') = Asifi^). (26) 
Po 

The RSch (/?2) dependence is obtained from the trun- 
cated Eqs. (I21D and (gT]) 



dAj{p^;f32) 1 d^Aj{p^;p2) 



dh 



2(31 d{\np?Y ' 



(27) 



where (j = 1, 2, . . .) and Ai = Ai- 

At loop- level four {n,^ = 4), where we include in RGE 
([T|) term with jmax = 5 (thus /^a), relations analogous to 
(EUl-dlll) can be found ^. 

It turns out that there is a clear hierarchy in magni- 
tudes |^i(g2)| > |^2(Q')| > |^3(Q')| > ••■ at all g^, 
in all or most of the anQCD models (cf. Fig. [5] for MA, 
Ml, M2; and Fig. [Din SecEUfor another model). 



loop-level = 4 , k,„^^ =6, A RSch 




FIG. 5: ^1 and A2 for various models (Ml, M2 and MA) with 
specific model parameters: cq = 2.94, Cr = 0.45, cj = 1.08 for Ml; 
Ci, = 0.1, Cp = 3.4 for M2; Uf = 3, A(„ 3) = 0.4 GeV in all three 
models. The upper three curves are Ai , the lower three are 3 X A2 ■ 
All couplings are in v-scheme (see Subsec. lVI Al l. A2 is constructed 
with our approach. Figure from: Ref. jlll |. 



We recall that the perturbation series of a space-like 
observable 'D{Q^) {Q^ = —q^ > 0) can be written as 



V{Q^ 



ipt 



a + dia + d2a + 



(28) 
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= fli + dia2 + \ d2- —di ) 03 H ,(29) 



where the second form (j29|) is the reorganization of the 
perturbative power expansion (|28p into a perturbation 
expansion in terms of a„'s (j24p (note: Oi = a). The 
basic analytization rule we adopt is the replacement 



An (n = l,2,...) 



(30) 



term- by-term in expansion (j29p . and this is equivalent 
to the analytization rule a" 1— > An term-by-term in 
expansion (|28p. However, in principle, other analyti- 
zation procedures could be adopted, e.g. a" 1— > ^", 
or a" I— > ^i,4„_i, etc. The described analytization 
a" I— > ^„ reduces to the MSSSh analytization in the case 

of the MA model (i.e., in the case of ^1 = ^j^ ), be- 
cause the aforementioned RGE-type relations hold also 
in the MA case. 

Let's denote by p("-)(g2) the TPS of ^ with 
terms up to (and including) the term ^ a"™, and by 

2?an° (Q^) the corresponding truncated analytic series 
(TAS) obtained from the previous one by the term- 
by-term analytization a" i— > An- The evolution of 
AkiQ'^) under the changes of the RSch was truncated 

in such a way that 92?!"" (Q^)/9/3j ~ -4n„+i (where 
j > 2). Further, our definition of ^^'s (fc > 2) via 
Eqs. (pi]) [cf. Eqs. (|^5|) ] involves truncated series which, 
however, still ensure the "correct" RScl-dependence 

SPan" (Q^)/9/i^ ^ ^n,„+i- This is all in close analogy 
with the pQCD results for TPS's: a2?("")(Q^)/5/3j ~ 

. In conjunction 



a"»+\ and dV^"'-'>{Q^)/dfi^ - a" 
with the mentioned hierarchy depicted in Fig. [51 this 
means that the evaluated TAS will have increasingly 
weaker RSch and RScl dependence when the number of 
TAS terms increases, at all values of Q^. 

On the other hand, if the analytization of powers were 
performed by another rule, for example, by the sim- 
ple rule a" ^^ Ai, the above RScl&RSch-dependence of 
the TAS would not be valid any more. An increasingly 
weaker RScl&RSch-dependence of TAS (when the num- 
ber of TAS terms is increased) would not be guaranteed 
any more. 



V. CALCULATION OF A^ FOR u NONINTEGER 

Analytization of noninteger powers in MA model was 
performed and used in Refs. [l5|, representing a gener- 
alization of results of Ref . 1231 . The approach was moti- 
vated by a previous work [26| where MA-type of analy- 
tization of expressions for hadronic observables was pos- 
tulated, these being integrals linear in a{tQ^) [similar to 
the dressed gluon approximation expressions, cf. Eq. (I44p 



and the first line of Eq. (|48)) ]. Analytization of noninte- 
ger powers a'^ or a'^ In a, is needed in calculations of pion 
electromagnetic form factor, and in some resummed ex- 
pressions for Green functions or observables, calculated 
within an anQCD model. 

In the mentioned approach, use is made of the Laplace 
transformation (/)l of fmiction / 



f{z) ^ if)Lit) : /(z) = / dte~^\f)Lit) , 

Jo 



where z = \n{Q^/A^). Using notations 
can be shown 



{an)L{t) = 
(An) Lit) = 



t" 



/3„"-^(n-l)! 

i3o-\n-iy. 



and (ESI), it 
(«)L(i) , (31) 



{Al)L{t) 



(32) 



Therefore, it is natural to define for any real v the fol- 
lowing Laplace transforms: 
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(A.) Lit) 
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(A)l(O • 



(33) 
(34) 



In MA model, at one-loop level, ia)Lit) and {Ai)Lit) are 
known 



a{z) = 
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Poz 



ia)Lit) = 
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M^^) - wl-^- 



{Ai)Lit) = j(^-T.Sit-k)j 



(35) 



(36) 



Since at one-loop A„ = A„, it follows in one- loop MA 
model 



A^z) 
Similarly, since 



dte-''——- 
/3or(j.) 



l-J2Sit~k) 



(37) 



fc=i 



a''(z)lna(z) = -—a'^iz) 
dv 



it can be defined 



dv 



a^iz) 



MA 



dv 



Ayiz) 



(38) 



To calculate higher (two-)loop level A^iz) in MA model, 
the authors of Refs. [ll| expressed the two-loop a(2)iz) 
in terms of one-loop powers aii\"^iz) In" an) (z) and then 
followed the above procedure. 



VI. EVALUATION METHODS FOR 
OBSERVABLES 

In pQCD, the most frequent method of evaluation of 
the leading-twist part of a space-like physical quantity 
is the evaluation of the available (RG-improved) trun- 
cated perturbation series (TPS) in powers of perturba- 
tive coupling a. Within the anQCD models, an analogous 
method is the aforementioned replacement a" i-^ An in 
the TPS (where An are constructed in Sec. lIV[) . and the 
evaluation thereof. More specifically, consider an observ- 
able T>{Q'^) depending on a single space-like physical scale 
Q^{= —q^) > 0. Its usual perturbation series has the 
form (P5)l . where a = a(^^; /32, /Sa, ■ • ■), with /x^ ~ Q^. For 
each TPS 2?(g2)|,'^^ of order TV, in the minimal anQCD 
(MA) model, the authors MSSSh [H, [13, \^ introduced 
the aforementioned replacement a" i— > A' 



(MA), 



p(Q2)W(MSSSh) ^ ^,(MA)^^^^^(MA)^. . . d^_^^^(MA) 

(39) 
This method of evaluation (via a" -^ An) was extended 
to any anQCD model in [13, [HI (cf. Sec.|W|. Further, in 
the case of inclusive space- like observables, the evaluation 
was extended to the resummation of the large- /3o terms: 



A. Large-/3o-niotivated expansion of observables 

We summarize the presentation of Ref. [ll| . We work 
in the RSch's where each (3k {k > 2) is a polynomial in 
Uf of order fc; in other words, it is a polynomial in (3q: 



k 



2,3,. 



(40) 



The MS belongs to this class of schemes. In such schemes, 
the coefficients dn of expansion (l28l) have the following 
specific form in terms of /?o : 

^(Q^)pt = a+(cii/3o + cio)a^ 

+ (c22/3o + C2i/3o + C20 + C2^-iPo^)a^ + ■■■ ■ (41) 

We can construct a separation of this series into a sum of 
two RScl-independent terms - the leading- /3o (L/3o), and 
beyond-the-leading-/?o (BL/3o) 



ppt = v: 



(L/3o) 
pt 



V 



(BL/3o) 
pt 



(42) 



where 



■,(L/3o) 
^pt 



Vir"" ^a + a^ [Pocii] + a^ [PoC22 + l3iCi 



^0^33 + ^/3o/3lC22 +/32C11 



IJ 



+ 0(/3o^a^). (43) 



Expression (j43|) is not the standard leading-/3o contribu- 
tion, since it contains also terms with Pj {j > 1), but only 
in a minimal way to ensure that the expression contains 



all the leading-/3o terms and at the same time remains 
RScl-independent. It can be shown that, for inclusive 
observables, all the coefficients in this L/3o contribution 
can be obtained, and can be expressed in the integral 
form [13 



p(W(Q2)^^^ r^F|(t)a(teCQ2) 
Jo '^ 



(44) 



where Fp(t) is the (Euclidean) L/3q -characteristic func- 
tion. In MS scheme, A = A which corresponds here to 
C = C = -5/3. No RScl /x^ appears in ([44]). Expression 
P^ is referred to in the literature sometimes as dressed 
gluon approximation. 

The BL/3o contribution is usually known only to ^^ a^ 



For it, we can use an arbitrary RScl fj. 



Q e ^ Q . Further, the powers a can be reexpressed 
in terms of an{^i^ 



a^ =a2-(/3i//3o)a3-f- 



flg- 



(45) 



Therefore, 



2?(Q')(TPS)=^^^^'''(Q')pt 



ht2 d2{Q^e^) + H asiQ'^e^) + U MQ^e^), (46) 

where ^2 = Cio is scheme-independent, and coefficients 
^3 and ^4 have a scheme dependence (depend on P2, Pi 
- i.e., on 62J and h^j). We note that expression ([46]) 
is not really a pure TPS, because its L/?o contribution 
(|43l) is not truncated. An observable-dependent scheme 
(D-scheme) can be chosen such that t^ — t^ — 0. For 
the Adler function V = d^, such a scheme will be called 
v-scheme. The analytization of the obtained 2?((5^)(tps) 
(|^S)) is performed by the substitution a„ 1— > An, Eq. (I30p . 
leading to the truncated analytic series (TAS) 



V{Q') = P(g2)(TAS) + 0(/3o'A) , (47) 

)(TAS) 

4«ioi'2(Q'e^')+i3^3(Q'e^Vi4A(Q'e^'). (48) 



p(g'^ 



^FUt)A,ite^Q') 
'0 ^ 



In the D-scheme, the last two terms disappear. Eq. ([48]) 
is a method that one can use to evaluate any inclusive 
space-like QCD observable in any anQCD model. As 
argued in Sec. IIVI the scale and scheme dependence of 
the TAS is very suppressed 



dV{Q^] 



(TAS) 



dX 



' PIA5 



'f3iA5 iX^lnfi^,Pj) . (49) 



If the BL/3o perturbative contribution is known exactly 
only up to (and including) ^a'^, then no t4 term appears 
in Eq. ([JS]) and the precision in Eqs. ([IT]) and ^^ is 
diminished: OiP^A^) >-^ 0{f3^Ai). 

It is interesting to note that the Taylor expansion of 
Aiite^Q'^) in V^■^'^o^Q'^)l,n in ([15]) around a chosen RScl 



TABLE I: Various order contributions to observables within 
PT, and MSSSh (=APT) methods ^,^: 



Process 


Method 


1st order 


2nd 


3rd 


GLS 
(Q - 1.76GeV) 


PT 
APT 


65.1% 
75.7% 


24.4% 
20.7% 


10.5% 
3.6% 


(M, = 1.78GeV) 


PT 

APT 


54.7% 
87.9% 


29.5% 
11.0% 


15.8% 
1.1% 



ln(/i^) reveals just the aforementioned a" i— > An analyti- 
zation of the large-/3o part (|43| . in any anQCD: 

= A1+A2 [PoCll] + A3 [/?oC22+/?lCll] 



+A4 



/5oC33 + 2 /3o/3lC22 +/32C11 



0(/3^A), 



where Ak = »4fc(^^;/?2,/?3, • ■ •)• I^^ other words, at the 
leading-/3o level, the natural analytization a 1— > ^1 in 
integral (|44p is equivalent to the terni-by-term analyti- 
zation a" !—)■ An [^ OLn ^-^■ ^n) in the corresponding 
perturbation series. This thus represents yet another 
motivation for the analytization a" 1— > An [<=> Eq. ([30]) 
postulated in Sec. |TV] of all the available perturbation 
terms in T). For the first motivation, based on the sys- 
tematic weakening of the RScl&RSch dependence of the 
truncated analytized P, see the end of Sec. IIVI 



B. Applications in phenomenology 

Evaluations in MA model, with the MSSSh-approach 
a" ^ A^^^ [H, [H, [3, are usually performed in MS 
scheme. The only free parameter is A (= A). Fitting the 
experimental data for T-decay, Z — > hadrons, e^e~ -^ 
hadrons, to the MSSSh approach for MA at the two- 
or three-loop level, they obtained A„^=5 « 0.26-0.30 
GeV, corresponding to: A„j.=3 w 0.40-0.44 GeV, and 
7r^^(MA)(-^2) j^ Q ^^24, which is above the pQCD world- 
average value as(M|) « 0.119 ± 0.001. The apparent 
convergence of the MSSSh nonpower truncated series is 
also remarkable - see Table ID 

In Refs. [10, llll, the aformentioned TAS evaluation 
method ^ in anQCD models MA (gl), Ml ^ and M2 
pop was applied to the inclusive observables Bjorken po- 
larized sum rule (BjPSR) db{Q'^), Adler function dv{Q'^) 
and semihadronic r decay ratio Vr The exact values of 
coefficients di and d? are known for space-like observ- 
ables BjPSR db{Q^) [2^ and (massless) Adler function 
dv[Q^) [2^, I30j. (The exact coefficient d^ of du has been 
recently obtained [3l| , but was not included in the analy- 
sis of Ref. [11] that we present here; rather, an estimated 
value of ^3 was used.) In the v-scheme, the evaluated 



TABLE II: Results of evaluation of rr{AS = 0,m, = 0) and 
of BjPSR di,(Q^) (Q^ = 2 and IGeV^), in various anQCD 
models, using TAS method (|48p . The experimental values 
are rr{AS = 0,mg = 0) = 0.204 ±0.005, db{Q'^ = 2 GeV^) = 
0.16 ± 0.11 and dt{Q^ = 1 GeV^) = 0.17 ± 0.07. 



dtjQ' = 2) 



dbjQ' = 1) 



MA 

Ml 

M2 



0.141 
0.204 
0.204 



0.137 
0.160 
0.189 



0.155 
0.170 
0.219 



massless dv{Q^) is 

dv{Q )(TAS) = 



dt 



^F,^(t)^i(te'^Q2;/32W,/33W) 



+ Y^A2{e^Q') 



(50) 



while BjPSR d6(Q )(tas) has one more term tsAsie^Q ). 
The difference between the (massless) true dx{Q^) (x — 

v,b) and dx{Q^){TAS) is 0{PqA4)- The semihadronic r 
decay ratio Vr is, on the other hand, a time-like quantity, 
but can be expressed as a contour integral involving the 
Adler function d^: 



rriAS^O,mq = 0) = 
2 /•"-' ds 
'0 



n 
1 

2^ 



rrir 



1 



Imn(s) 



+ 71- 



(1 + e"^f{l - e'^)d^{Q^ - rur^e'^). (51) 



This implies for the leading-/3o term of r,- 



r^(AS' = 0,m, = 0) 



_m(L/3o) _ 



jF,^it)^,ite^m^% 



(52) 
where 2ti is the time-like coupling appearing in Eqs. ([5])- 
([5]), and superscript A4 in the characteristic function in- 
dicates that it is Minkowskian (time-like). The latter was 
obtained by Neubert (second entry of Refs. [23] )■ The 
beyond-the- leading- /3o (BL/3o) contribution is the con- 
tour integral 



(BL/3o) ^ 

(l + e»'^)3(i_e*'^)l2(e^m,2e*^). (53) 



r^(AS' = 0,mq = 0) 
1 
24^ 

The parameters of anQCD models Ml ^ and M2 HH) 
were then determined [lH by fitting the evaluated ob- 
servables to the experimental central values rT-(A5 = 
0,m, = 0) = 0.204 (for Ml and M2), and to db{Q^ = 
IGeV^) = 0.17 and 4(Q^ = 2) = 0.16 (for Ml). For Ml 
we obtained: c/ = 1.08, Cr = 0.45, cq = 2.94. For M2 we 
obtained: c„ =0.1 and Cp — 3.4. 

The numerical results were then obtained [11| . In mod- 
els MA, Ml and M2 they are given for r^ in Table HII for 




- 

-0.1 
-0.2 



loop-le\'el=4, k^^^B 

A= 0.4 GeV 

M1:Co=2.94, c,=0.45, C(=1.08 

M2:(^=0.1,Cp=3-4 

f MA,M1,M2: ji-=Q^e"^^(masslessp.] 
mass corr. and ST PS; fi'=0^ 



STPS4-l.(bWS) 

STPS3-l.(bMS) 

M2{A) 

Ml (A) 

MA (A) 

exp. 



0.4 



0.6 



1 



1.4 



1.6 



1.1 



1.2 
Q [GeV] 
FIG. 6: Adler function as predicted by pQCD, and by our ap- 
proach in several anQCD models: MA, Ml, M2. The full quantity 
is depicted, with the contribution of massive quarks included. The 
experimental values are from |32l |. Figure from: Ref. [lll |. 



Adler function ^^(Q^) in Fig. El and for BjPSR ^(Q^) 
(in Ml and M2) in Figs. [7] and M (Table M and Figs. [3 
[3 [8] are taken from Ref. [11|). All results were calculated 
in the v-scheme. For details, we refer to Ref. [11[|. 



0.17 
0.16 
0.15 

^ 0.13 

0.12 

0.11 

0.1 




|M|I 

M1:C|j4l94, c,-b.4d c,=1 



1 



1.5 



3.5 



2 2.5 3 

Q [GeV] 

FIG. 7: Bjorkcn polarized sum rule (BjPSR) di,{Q^) in model 
Ml, in various RSch's and at various RScl's. The vertical lines 
represent experimental data, with errorbars in general covering the 
entire depicted range of values. 

Analytic QCD models have been used also in the 
physics of mesons [33|, |3J| , in calculating various meson 
masses by summing two contributions: that of the con- 
fining part and that of the (one-loop) perturbative part 
of the Bethe-Salpeter potential. In Refs. [33], the (one- 
loop) MA couplin g [3] was used to calculate/predict the 
masses; in Refs. [3J|, the experimental mass spectrum 
was used to extract the approximate values of the (ana- 
lytic) coupling Ai{Q^) at low Q^. In this formalism, the 
current quark masses were replaced by the constituent 
quark masses, accounting in this way approximately for 
the quark self-energy effects. The results by the authors 
of Ref . ,34| indicate that Ai{Q'^) remains finite (and be- 
comes possibly zero) when Q^ — > 0. 



O 



0.2 

0.19 
0.18 
0.17 
0.16 
0.15 
0.14 
0.13 
0.12 



-T^-^ 



(a) ti 



loop-level-3, \^-5 
model M2; c„^0.1 Crt-3 4 



LS(A) 

LS(B) 

n'=0^ (A) 

li'=Q^ (B) 

(A) 

(B) 

sk(A) 



M2 



A= 0.4 SeV 



1 



1.5 



3.5 



2 2.5 3 

Q [GeV] 

FIG. 8: As in the Fig.[3 but this time for model M2. Both figures 
from: Ref. [ull. 



VII. ANALYTIC QCD AND ITEP-OPE 
PHILOSOPHY 

In general, the deviations of analytic ^i(Q^) from 
the perturbative coupling apt{Q^) at high Q^ ^ A^ are 
power terms 

k 



4) (Q^»A^), 



\SAiiQ^)\ = \Ai{Q^)-apt{Q^ 



where fc is a given positive integer. Such a coupling intro- 
duces in the evaluation (of the leading- twist) of inclusive 
space-like observables V{Q'^), already at the leading- /3o 
level, an UV contribution SV^^^^Q^) which behaves like 
a power term [l8| 



5I?(UV)(g2) 



A 2 \ min(k,n) 



[{ k ^ n , 



(54) 



where n e Af is the position of the leading IR renormalon 
of the observable V{Q^): ii k = n, then the left-hand side 
of Eq. dMD changes to (AVQ^)" ln(AVg2) [l|. Such 
nonperturbative contributions coming from the UV sec- 
tor contradict the ITEP Operator Product Expansion 
(OPE) philosophy (the latter saying that such terms can 
come only from the IR sector) |35| . 

Two specific sets of models of anQCD have been in- 
troduced in the literature so far such that they do not 
contradict the ITEP-OPE: 

(A) a model set based on a modification of the (3{a) 
function jl7|: 

(B) a model set obtained by a direct construction [l^ . 



A. Set of models A 

This is the set of models constructed in Refs. iTTj]. The 
TPS /3(a) used in pQCD is 



da 



N 



51nQ2 



/3W(a) = _/?oa2 | l + J^Cja^ 
j=i 



(55) 



This was then modified, /?'^^^(a) i-^ /3'^^^(a), by fulfiUing 
three main conditions: 

1.) P^^' (a) has the same expansion in powers of a as 

2.) /3(^)(a) CaP with C > and p < 1, for a > 1, 

in order to ensure the absence of Landau singularities; 

3.) /3*-''^''(a) is analytic function at a = 0, in order to 
ensure \a{Q^) - apt{Q^)\ < (A^/Q^)'" for any fc > at 
large Q^ (thus respecting the ITEP-OPE approach). 

This modification was performed by the substitution 
a |—^ u{a) = 0/(1 + 770), r] > being a parameter, and 



/3(^)(a) = -po 



N 



K{a~u{a)) + y^ Cju{a) 



J+2 



3=0 



,(56) 



and Cj are adjusted so that the first condition is fulfilled 
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FIG. 9: The couplings A2 and ^3, together with the corresponding 
coupling yli , are plotted as a function of Q, in the MS-scheme, with 
A = 0.4 GeV. The parameters used for the couplings are -r) = 0.3, 
hi = 0.1, and ^2 = 0. Figure from: Ref. [l|. 



Co = 1-??K, 



Cl 



Ci + 2ri~r/ K, etc. 



This procedure results in an analytic coupling a{Q^), 
withp = 1 and ( = PqK, and with two positive adjustable 
parameters k and ry. The QCD parameter A was taken 
the same as in the pQCD. Evaluation of observables was 
carried out in terms of power expansion, with the replace- 
ment ap\ I— > a". Further, the couplings in this set are IR 
infinite: a{Q'^) - l/iQ'^f^'' -^ 00 when Q^ ^ q. These 
new a((5^)'s are analytic (a = ^1). The RScl and RSch 
sensitivity of the modified TPS's of space-like observables 
turned out to be reduced. The author of Refs. [l3| chose 
K = l//3o; by fitting the predicted values of the static in- 
terquark potential to lattice results, he obtained 77 « 4.1. 



B. Set of models B 



in the UV region, the role of the perturbative coupling. 
Lo and Li are analytic and chosen aiming at a low km- 
dependence in the IR region. 



1 



1 

T 



^i^i(l-V^) 



1 



-5,(2;), y^>Q, i = 0, 1. (58) 



Functions gi{x) are chosen in simple meromorphic form 
2x 



9oix) 
giix) 



(l + i/o)+a;(l-t'o)' 



d + X 



< j/o < 1; (59) 
rf > 0, (60) 



with the constants fixed at typical values i^q = 1/2 and 
i^i ~ d = 2. The additional expoinential term in (|57p is 
chosen in a similar meromorphic form 



This is the set of models for Ai constructed in Ref. [T^| . 
A class of IR-finite analytic couplings which respect 
the ITEP-OPE philosophy can be constructed directly. 
The proposed class of couplings has three parameters 
{•q,hi,h2). In the intermediate energy region (Q ~ 1 
GeV), the proposed coupling has low loop- level and 
renormalization scheme dependence. We outline here 
the construction. We recall expansion ^ for the per- 
turbative coupling a{Q'^), where L = logQ^/A^ and Km 
are functions of the /3-function coefficients. This expan- 
sion (sum) is in practice usually truncated in the index 
k (k < km)- The proposed coupling is obtained by modi- 
fying (the nonanalytic) i's to analytic quantities Lq and 
Li that fall faster than any inverse power of Q^ at large 
Q^, and by adding to the truncated sum another quantity 
with such properties: 






fe=i e=o 



Tk 
^0 



"v^/(x), (57) 



where x = Q^/A^. The second term is only relevant 
in the IR region, and the first term (double sum) plays. 



r"^f{x)^h. 



1 



(H-a;/2) 



l^f,-vV^ 



(61) 



Results for Ai, A2 and A3, for specific typical values of 
parameters 77, hi and /12, are shown in Fig. [51 Couplings 
A2 and ^3 are constructed via A2 and .43, according to 
the procedure described in Sec. IIVI Eqs. (P5)l . 

A general remark: if Ai{Q'^) differs from the pertur- 
bative a{Q^) by less than any negative power of Q^ at 
large Q^ (^ A^), then the same is true for the difference 
between any Ak{Q^) and ak{Q^) (fc = 2,3, . . .). 



VIII. SUMMARY 

Various analytic (anQCD) models, i.e., analytic cou- 
plings Ai{Q^), were reviewed, including some of those 
beyond the minimal analytization (MA) procedure. 

Analytization of the higher powers a" 1-^ An was con- 
sidered; an RGE-motivated approach, which is applicable 
to any model of analytic Ai , was described. Analytiza- 
tion of noninteger powers a'^ in MA model was outlined. 
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Evaluation methods for space-like and time-like ob- 
servables in anQCD models were reviewed. A large-/?o- 
motivated expansion of space- like inclusive observables is 
proposed, with the resummed leading-/3o part; on its ba- 
sis, an evaluation of such observables in anQCD models is 
proposed: truncated analytic series (TAS). Several eval- 
uated observables in various anQCD models were com- 
pared to the experimental data. We recall that evaluated 
expressions for space-like observables in anQCD respect 
the physical analyticity requirement even at low energy, 
in contrast to those in perturbative QCD (pQCD). 

Finally, specific classes of analytic couplings Ai{Q'^) 
which preserve the OPE-ITEP philosophy were dis- 



cussed, i.e., at high Q^ they approach the pQCD cou- 
pling faster than any inverse power of Q^. Such analytic 
couplings should eventually enable us to use the OPE 
approach in anQCD models. 
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